Tight-binding parameters for charge transfer along DNA 
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We systematically examine all the tight-binding parameters pertinent to charge transfer along 
DNA. The tt molecular structure of the four DNA bases (adenine, thymine, cytosine, and gua- 
nine) is investigated by using the linear combination of atomic orbitals method with a recently 
introduced parametrization. The HOMO and LUMO wavefunctions and energies of DNA bases 
are discussed and then used for calculating the corresponding wavefunctions of the two B-DNA 
base-pairs (adenine-thymine and guanine-cytosine) . The obtained HOMO and LUMO energies of 
the bases are in good agreement with available experimental values. Our results are then used for 
estimating the complete set of charge transfer parameters between neighboring bases and also be- 
tween successive base-pairs, considering all possible combinations between them, for both electrons 
and holes. The calculated microscopic quantities can be used in mesoscopic theoretical models of 
electron or hole transfer along the DNA double helix, as they provide the necessary parameters for a 
tight-binding phenomenological description based on the tt molecular overlap. We find that usually 
the hopping parameters for holes are higher in magnitude compared to the ones for electrons, which 
probably indicates that hole transport along DNA is more favorable than electron transport. Our 
findings are also compared with existing calculations from first principles. 

PACS numbers: 73.63.-b, 82.39. Jn,87.14.gk, 87.15.A- 



I. INTRODUCTION 

DNA plays a fundamental role in genetics and molecu- 
lar biology since its sequence of bases, adenine (A), gua- 
nine (G), cytosine (C), and thymine (T), contains the 
genetic code of living organisms. During the last decade, 
DNA and its charge transport properties have attracted 
the interest of a large interdisciplinary community due 
to its potential use for nanodeviceSjeither for assembling 
nanocircuits or as a molecular wire [J 0, d, 01 • Charge mi- 
gration through DNA could also play an important role 
in biology; it may be a critical issue in carcinogenesis and 
mutagenesis 0, |g] . For example the rapid hole migration 
from other bases to guanine is connected to the fact that 
direct strand breaks occur preferentially at guanines 
Long-range charge transfer along the 7r-stacking of the 
DNA double helix may also be crucial for DNA damage 
and repair 0, Q . 

However, it is not yet clear "how much" DNA con- 
ducts. Experiments cover a wide range of behavior. In 
particular, previous results found that A-DNA covalently 
bonded to Au electrodes is an insulator Q. Furthermore, 
insulating behavior -both for single DNA molecules as 
well as for small bundles of DNA molecules- was observed 
at the 100 nm length scale this was confirmed both 
for mixed base-pair sequence as well as for homogeneous 
poly(dG)-poly(dC), for lengths between contacts in the 
range 40-500 nm, for substrate Si02 or mica, and for elec- 
trode material gold or platinum. The absence of dc con- 
ductivity in A-DNA was also reported for DNA molecules 
adsorbed on mica, in agreement with the authors' first 
principles electronic structure calculations (Tl| . 

In a different direction point the results of 
Yoo et al. [12], who investigated the electrical transport 
through poly(dA)-poly(dT) and poly(dG)-poly(dC) 



DNA, i.e. molecules containing identical base-pairs. 
Their experimental results suggest that electrical trans- 
port through DNA molecules occurs by polaron hopping 
and the possibility of a DNA field-effect transistor 
operating at room temperature was also demonstrated 
Moreover, gate- voltage dependent transport mea- 
surements showed that poly(dA)-poly(dT) behaves as 
a n-type semiconductor, whereas poly(dG)-poly(dC) 
behaves as a p-type semiconductor [12l |. In another 
report it has been shown that poly(dG)-poly(dC) can 
act as a semiconducting nanowire exhibiting better 
conductance than poly(dA)-poly(dT) [13]. 

Measurements of electrical current as a function of the 
potential applied across single DNA ropes at least 600 nm 
long, indicated metallic-like gapless behavior and efficient 
conduction comparable to that of conducting polymers 
[l^ . In this work it was mentioned that the observed 
behavior together with the fact that DNA molecules of 
specific composition and length, ranging from a few nu- 
cleotides to several tens of /im, can be prepared, makes 
DNA ideally suited for the construction of mesoscopic 
electronic devices [ijj. 

The ac conductivity of DNA has been also explored. 
For example, A-phage DNA at microwave frequencies 
showed strong temperature dependence of the ac conduc- 
tivity around room temperature, exhibiting a crossover 
to a weak temperature dependence at lower tempera- 
tures [IB]- Measurements of the quasi-static, frequency- 
dependent conductivity below 1 MHz performed on wet- 
spun, macroscopically oriented, calf thymus DNA bulk 
samples, showed that the electrical conductivity can be 
rather well described by an activated Arrhenius law with 
activation energy of 0.9 eV or by hole hopping p^ . 

The observed important deviations in the conductivity 
of DNA can be -at least partially- attributed to several 
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external conditions like the type of the substrate, the dis- 
tance between the electrodes, and the contact material, 
in conjunction with the soft nature of the biomolecule. 
In addition, several intrinsic characteristics, like the lo- 
cal chemical environment which is a result of the solu- 
tion involved in the DNA preparation the hydrogen 
bonding fis'l, and the degree of stretching of DNA |19j 
affect the electronic properties. 

The great effort to ameliorate DNA's transport prop- 
erties led to alternative molecular conductive candidates 
based on DNA. In M-DNA the imino-proton of each base- 
pair has been substituted with a metal ion [23|. In con- 
trast to B-DNA, M-DNA presents no plateau in its I-V 
curve and exhibits metallic-like conduction at room tem- 
perature [1^. Size-expanded DNA bases are produced 
by the addition of a benzene ring to a natural DNA base 
[2l| and have been combined with natural bases to form 
xDNA and yDNA, a new class of synthetic nucleic acids. 
It has been theoretically shown that xDNA and yDNA 
have smaller HOMO-LUMO gap than B-DNA [Ij, i.e., 
they might function as molecular wires better than nat- 
ural DNA. 

Many mechanisms have been suggested to explain the 
charge transport along DNA, where generally 7r-pathway 
transfer due to the overlap of tt molecular orbitals of the 
stacked aromatic bases of DNA can lead to charge prop- 
agation even at long distances i, 0, Hi. Examples of 
mechanisms that have been examined are band trans- 
port [H, , polaronic transport [H, [25I .J 
fluctuation-facilitated charge migration [30|, |3l|, |32, 
and variable range hopping [s^]. 

The aim of this work is to provide electronic param- 
eters for charge (electron or hole) transfer along DNA, 
assumed that the transport mechanism is the 7r-pathway. 
Using the linear combination of atomic orbitals (LCAO) 
method, our investigation focuses on the calculation of 
HOMO (highest occupied molecular orbital) and LUMO 
(lowest unoccupied molecular orbital) vr molecular wave- 
functions and energies for the four DNA bases. A, G, C, 
T, and the two B-DNA base-pairs, adenine-thymine (A- 
T) and guanine-cytosine (G-C). Then the hopping matrix 
elements between neighboring bases or base-pairs, for all 
possible combinations between them, are provided, for 
both electrons and holes. The obtained hopping param- 
eters and the HOMO and LUMO energies can be used 
in a tight-binding phenomenological description of charge 
transfer along DNA. Such a mesoscopic description is ap- 
plied in numerous theoretical models of charge transport 
In m, [13, H, m H m, m, [13] . our results are com- 
pared with existing calculations from first principles, as 
well as with available experimental data for the HOMO 
and LUMO energies of DNA bases. 

The article is organized in the following way: In Sec. [IT] 
we explain in detail our theoretical approach. Specifi- 
cally, in Sec. Ill Al we apply a simple LCAO method for 
the TT electronic structure of the four DNA bases using 
a novel parametrization [s^, in Sec. Ill Bl we discuss the 
HOMO and the LUMO of base-pairs by employing a sim- 



ilar linear combination of molecular orbitals method, and 
finally Sec. Ill CI is devoted to the tight-binding descrip- 
tion and the calculation of charge transfer parameters in 
this framework. In Sec. IIIII we present and discuss our 
results, along with a comparison with earlier theoretical 
calculations as well as with available experimental values. 
Finally, in Sec. IIVI we state our conclusions. 



II. THEORY 

A. TT molecular structure of DNA bases 

In this subsection we present the LCAO methodol- 
ogy used for the description of the tt electronic struc- 
ture of the four isolated DNA bases. A, G, C, and T. 
DNA bases are planar organic molecules bonded by sp^ 
hybridization, where the atoms have their pz atomic or- 
bitals perpendicular to the molecular plane. The elec- 
trons that occupy these atomic orbitals are delocalized 
to form TT molecular orbitals. The LCAO method pro- 
vides a straightforward approach to obtain the tt molec- 
ular structure. In its simpler form, used in the present 
article, a tt molecular single-electron wavefunction can be 
approximated as 



N 



^\r) = J2c.pl{r) 



(1) 



The index i implies summation over all atoms (N totally) 
which contribute pz electrons to the DNA base of inter- 
est. |cip gives the probability of finding the electron oc- 
cupying the molecular orbital \E'^(r) at the atom indexed 
by i, while pl{r) denotes the corresponding atomic, or 
atomic-like, orbital. The molecular wavefunction obeys 
Schrodinger equation H'''ii''{r) = £'''\E'^(r), where E'' is 
the eigenergy of the base. Substituting in the latter equa- 
tion ^''(r) through Eq. ([1]), multiplication with the con- 
jugate orbital *(r), and integration over all space, gives 
a linear system of A'^ equations obeyed by the unknown 
coefficients Ci of the molecular wavefunction under deter- 
mination and its energy eigenvalue E^. Following a stan- 
dard procedure and assuming that the atomic or atomic- 
like Pz orbitals located at different atoms are orthogo- 
nal, one obtains that the solution of the linear system 
of equations determining Ci and the corresponding E'^, 
is equivalent to the diagonalization of the Hamiltonian 
with matrix elements H^j = J d^r pl*{r) pH^) (see 
for example Ref. [ssj). 

Through numerical diagonalization of the symmetric 

N X N Hamiltonian matrix i/f!, we obtain the coefficients 

''J 

Ci (eigenvectors) providing the molecular orbitals ^'''(r) 
[cf. Eq. ID)], as well as the corresponding eigenenergies 
E^. To this end we need the values of the Hamiltonian 
matrix elements, iJj'j. We employ a recently introduced 
parametrization that has been successfully used for the 
description of the energy of tt frontier orbitals in more 
than sixty planar organic molecules [ssl], including the 
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DNA bases, as well as in flavin [39[, and has shown 
to provide accurate calculations of the ttHOMO ener- 
gies and TT-TT* gaps. Regarding the diagonal matrix ele- 
ments, -fffj = Ei, we use Ec — —6.7 eV for carbon atoms, 
EN2 = —7.9 eV for nitrogen atoms contributing one Pz 
electron (i.e. with coordination number 2), en-^ — —10.9 
eV for nitrogen atoms with two pz electrons (i.e. with 
coordination number 3), and eq — —11.8 eV for oxygen 



(2) forming the base-pair. Inserting Eq. ([3]) into 
Schrodinger equation H''p *^/^(r) = e''^^^ "^h/tM) 



'H/L ^ H/L^ 



(where e'^^^^ denotes the HOMO/LUMO base- 
pair energy), multiplying once with ^'^J^*(r) 
and once with 5'^^^*(r), and integrating over 
ah space, using ^'^H(r) = T,^=icf(uP'z^\r), 



atoms. For the nondiagonal matrix elements H^j (i / j) \lr^^J' (r) — J2^^i c^/^p'z'^^ (r) [cf Eq ([T])] we ob 



referring to neighboring sp^-bonded atoms we use the 
following expression proposed by Harrison : 



ppn 



-0.63- 



(2) 



where m is the electron mass and d is the distance be- 
tween the corresponding nearest-neighboring atoms. All 
other nondiagonal matrix elements, referring to atoms 
not sp^-bonded are assumed equal to zero, H^j = 0. 

Diagonalizing the Hamiltonian matrix, one obtains N 
molecular orbitals and their eigenenergies. The lower en- 
ergy orbitals are successively filled by two electrons each, 
until to accommodate all available electrons. Then the 
highest-energy occupied orbital, \E'^(r), is the tt HOMO 
and the lowest-energy unoccupied one, 5'|^(r), is the tt 
LUMO. Following the usual convention, we denote oc- 
cupied molecular orbitals by tt and unoccupied ones by 



tain the following system of equations 



pb(i) r 



- i^H/L '-2 



E 



E 



hp 

H/L 
hp 

H/L 



Cl 
C2 



(4) 



Here we have assumed that Pz orbitals belonging to 
atoms of different bases are orthogonal and also that 

J H/L 



E 



b{m) 



H/L 



fj^j^, for m=l or 2. The overlap integral, tu/L = 
/d^r -^ffj* H'P ^f/^, which expresses the hopping 
parameter for a charge (hole/electron) transfer between 
the two bases of the base-pair, equals 



EE 

i=i j=i 



H/L * H/L y 
S(l) S(2) ^^3' 



(5) 



B. HOMO and LUMO of B-DNA base-pairs 

Regarding the B-DNA base-pairs, we follow a differ- 
ent procedure in order to obtain the HOMO and LUMO 
wavefunctions. This is because the two bases (A and T, 
or G and C) are connected with non-covalent hydrogen 
bonds to form the base-pair (A-T or G-C, respectively). 
The length of hydrogen bonds, around 3 A, is longer than 
the typical length of a covalent bond connecting neigh- 
boring atoms within a base, which is around 1.3-1.5 A. 
In other words, a base-pair is not assumed in our calcu- 
lation like a single molecule but, instead, it is treated like 
two adjacent molecules with electronic overlap. However, 
we still use the terms HOMO and LUMO for the base- 
pairs, meaning the single-electron wavefunctions which 
represent the highest in energy occupied orbital and the 
lowest in energy unoccupied orbital, respectively, of the 
molecular complex. We assume that these wavefunctions 
describe an inserted hole or electron, respectively, within 
a base-pair. Below we follow a straightforward linear 
combination of molecular orbitals approach to designate 
the base-pairs' HOMO and LUMO. Namely, the base- 
pair HOMO/LUMO {H/L) wavefunction reads 



H/L 



m 
H/Ly 



(3) 



,6(2 



where 5'^^^(r), ^'^^^(r) are the corresponding 



- H/L 



where = f d^r pl^^^ *(r) H^p pi'^^\r). 

Because the base wavefunctions \['''(r) are real, the co- 
efficients Ci in Eq. ([T]) are real, and the same holds for 
the overlap integrals tn/L, i-e., ^jj/^ =tH/L- The matrix 
elements Vij are generally provided through the Slater- 
Koster expression (4l| : 



(6) 



where (j) denotes the angle formed by the line connecting 
atoms i and j and the plane perpendicular to Pz orbitals 
(i.e. the plane of bases). More details regarding the 
meanin g of Vppa and Vpp.^ matrix elements can be found 
in Ref. li^. 

For atoms belonging to different bases within a base- 
pair, i.e. the case of interest in this subsection, the angle 
(f) — Q and Vij = T^p^r • Note that for intra-base covalently 
bonded neighboring atoms, Ipp^r ^-^^ the interatomic ma- 
trix elements proposed by Harrison [iQl, cf. Eq. 
However, Harrison's relations are valid for interatomic 
distances of the order of those of covalent bonds. When 
dealing with larger interatomic distances, e.g. between 
atoms belonging to different molecules, the proportional 
to 1/d^ expressions of Harrison [l^l are replaced by ap- 
propriate exponentially decaying expressions of the form 



* ppiT — : 



(7) 



HOMO/LUMO orbitals of the two bases (1) and 



where the constants A, (3 are determined through the 
requirement that at a typical covalent bond distance, d^, 
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the values of this expression and its derivative in respect 
to d, coincide with those of Harrison's expressions given 
in Eq. This yields A = -OMfi^/mdl and f3 = 2/do- 
Here we choose do = 1-35 A. 

Using Eq.© for Vij (since Vij is equal to Ipp^ for 
atoms of different bases within the same base-pair, i.e. 
(j) ~ 0) and the known coefficients of the base 

HOMO/LUMO (obtained as described in the previous 
subsection), then tfj/L are calculated through Eq. (O. 

The quantities and E'^j'^f^ in Eq. (g]) are the 

HOMO/LUMO eigenencrgies of the corresponding bases. 
Therefore, the 2x2 system of Eq. (|3|) can be analytically 
solved to determine Ci, C2, and E^^j^. To obtain the 
HOMO (LUMO) of the base-pair, the higher (lower) en- 
ergy solution of the 2x2 system Q is considered and E^^ 
{e'^) is the corresponding eigenenergy. From the values 
Ci and C2 of the solution, the base-pair HOMO /LUMO 
wavefunction is obtained through Eq. ([3]). For later use 
it is mentioned that by multiplying Ci and C2 with the 
LCAO coefficients of the corresponding base wavefunc- 
tions ^"^/^ and ^'^^^ [cf. Eq. ([1])], the base-pair wave- 
function ^ can be equivalently written as 



N 



1=1 



(8) 



Here the sum is extended over the N atoms (contributing 
Pz electrons in tt bonds) of the whole base-pair {N = 18 
for A-T, while iV = 19 for G-C). 



Tight-binding parameters for charge transfer in 
B-DNA 



HOMO and LUMO energies of bases (or base-pairs) as 
well as hopping parameters between successive bases (or 
base-pairs) -calculated through a similar procedure like 
in the previous subsection, cf. Eq. ([5])- provide an esti- 
mate of the parameters used in tight-binding models that 
describe charge transport along DNA. Such phenomeno- 
logical models facilitate larger scale simulations, reaching 
a mesoscopic level of description. Depending on the par- 
ticular problem of interest, a tight-binding description 
at the base-pair level, or at the single-base level, may 
be more appropriate. The relevant parameters for both 
cases are discussed below and their estimated values are 
presented in Section HL 



1. Description at the base-pair level 

A tight-binding description of charge transfer between 
successive base-pairs A — 1, A, A -I- 1, ... of double- 
stranded DNA can be obtained considering that extra 
electrons inserted in DNA travel through LUMOs, while 
inserted holes travel through HOMOs. In this approxi- 
mation the time-dependent single carrier (hole/electron) 



wavefunction of the whole macromolecule, ^^'^j^{r, t), is 
considered as a linear combination of base-pair wavefunc- 
tions with time-dependent coefficients, i.e.. 



A 



(9) 



where *^/^^r) is the A*'' base-pair's HOMO/LUMO 
wavefunction and the sum is extended over all base-pairs 
of the DNA molecule under consideration. 

Starting from the time-dependent Schrodinger equa- 

tion, ih—^ = H^'^^-i'^fj^, and following a similar 
procedure and assumptions similar to the ones of the pre- 
vious subsection, one obtains that the time evolution of 
the coefficients ^a(^) obeys the tight-binding equations 



ih- 



dt 



E 



,6p(A) 
H/L 



. 6p(A:A-l). bp(A;A+l), 



(10) 

Here, iJ^^^^^ is the HOMO/LUMO energy of base-pair 
A, as obtained in the previous subsection IIIBI Depend- 
ing on the actual sequence of the DNA molecule under 
consideration, two values are possible for the on-site en- 
ergies E^lf^ of Eq. (Uni), corresponding to A-T or G-C 
base-pairs. Using the HOMO or LUMO base-pair wave- 
functions (^"^ or cf. Eq. ([5]), the corresponding 
hopping parameters (i^ or i'l^ between successive base- 
pairs are obtained by 



MA;A') 
HjL 



i=l j=l 



(11) 



where the indices A, A' denote neighboring base-pairs and 



The sums over i 



and j in Eq. pT|) extend up to the total number of atoms 
Nx and Nx> , respectively, constituting the corresponding 
base-pair [cf. Eq. ([5])]. Here there is a difference with the 
hopping integrals in Eq. ([5]), since in the latter case the 
sums up to Ni and N2 extend over the total number of 
atoms constituting the corresponding bases. The matrix 
elements Vij in Eq. (jlip are given by the semi-empirical 
Slater-Koster expression, Eq. ([5]). Now (j) ^ and T^p^ is 
obtained through Eq. ([7]), as before. Vppa is also obtained 
from the same expression, Eq. ([7]) but with a different 
preexponential coefficient, A — 2.22h? /md^, resulting 
from the constant 2.22 which appears in the correspond- 
ing to Eq. Harrison's formula for Vppa [i^l- There- 
fore. using standard expressions from solid state physics 
[40I I411 143 . 143I , we are able to calculate the interatomic 
matrix elements Vij . The Slater-Koster expression of Eq. 



® (see Table I of Ref. j4i]) with matrix elements V^ 



ppTT 



and Vppa of the Harrison's type [40|, is typically used in 
solid state physics and it has been also applied in the 
case of DNA charge transfer in Ref. 0]. Our methodol- 
ogy for calculating the overlap matrix elements is very 
similar to that of the latter work. The difference is that 
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we use the exponentially decaying expressions of the form 
of Eq. ^ for Ypp^^ and Vpp„, adopted from Refs. [4^.[43j. 
instead of the fitting of similar expressions with ab ini- 
tio calculations as used in Ref. AH coefficients and 
matrix elements in Eq. are real, resulting in real 

bp(A;A') _ ,bp(A';A) 
''H/L " ''H/L 

Knowledge of the coefficients C^^^ from the base-pair 
wavefunctions of Sec. Ill Bl [cf. Eq. ([5])] and of the ma- 
trix elements Vij from the geometrical structure of DNA 
and the formulae discussed above, allows the calculation 

of the hopping integrals t''^j^^ ^ [cf. Eq. pTjl ] which ap- 
pear as parameters in the tight-binding Eq. (fTU)) . There- 
fore, one can use the tight-binding parameters E'^^j^ and 

^ ' computed in this work (cf. Table [Til and Table [1111 
respectively, in Sec. lIIII) . in order to numerically solve the 
system of equations (ITU)) and obtain, through Ax{t), the 
time evolution of a charge [cf. Eq. ([9])] propagating along 
any DNA segment. 



ence is that the carrier wavefunction of the whole macro- 
molecule is now considered as a linear combination of 
single base wavefunctions, instead of base-pair wavefunc- 
tions. Therefore 



(12) 

where A denotes base-pairs, the sum is again over all 
successive base-pairs of DNA, and ^^h/'l 
HOMO/LUMO wavefunctions of bases in the A*'* base- 
pair, located at the one and the other DNA strands, re- 
spectively. 



2. Description at the single-base level 

The tight-binding description at the single-base level 
is similar to that in the previous subsection. The differ- 



Under the same assumptions as previously and consid- 
ering only neighboring hoppings (between adjacent base- 
pairs), the tight-binding equations for the time depen- 
dent coefficients of Eq. read 



ih 
ih 



dt 
dt 



E 



6(A,1) 
H/L 



A. 



6(A.1;A,2) „ 



6(A,1;A-14) 4 , - 



E 



H/L 



B 



A + 



6(A,2;A,1) 
H/L 



^> + ^H/L ^^-i 



H/L 

,b(A,2:A+l,2) „ 
^H/L -^^+1 



6(A,1;A-1,2) 
H/L 



,6(A,1;A+1,2) „ 
t^/L 



6(A,2;A-1,1) . b(A,2:A+l,l) . , „^ 

%/L ^^-1 + ^HIL AxilA) 



H/L 



Here E^^^^ are base HOMO/LUMO energies and the 
hopping parameters i^/^ are interbase transfer inte- 
grals of the general form of Eq. ([5]) with Vij given by 
Eq. ([5]). = for interbase transfer integrals within the 
same base-pair (interstrand intra-base-pair hoppings), 
^b(A^i,A,2) _ ^b{^\,2,\,i) ^ while (/) ^ for interbase trans- 
fer integrals between successive bases within the same 
strand (intrastrand hoppings), or between 

diagonally located, at different strands, bases (inter- 
strand inter-base-pair hoppings), t^^j^''''^^'"'^ with i ^ j. 
VppTT and Vpp„ are given through exponentially decaying 
Eq. ([7]), as discussed before. 

A tight-binding description at the single-base level re- 
quires solving of a double number of differential equa- 
tions, Eq. (fT^ . compared to the description at the base- 
pair level, Eq. p^ . However, it allows individual base 
properties to be taken into account when necessary. Sec- 
tion |TTT] contains the tight-binding parameters which ap- 
pear in Eq. ^] see Table |T] for E\j,^, Table [IV] for 



iffji Table|V]for and Tables |Vl] and |VlT] 

for t^^^''^^"^'"'^ with i ^ j. 

III. RESULTS AND DISCUSSION 

The results and their discussion have been organized in 
two subsections. In subsection lHI Al we present our calcu- 
lations for the HOMO and LUMO energies as well as for 
the corresponding wavefunctions of the four DNA bases, 
adenine, guanine, cytosine, and thymine, and of the two 
base-pairs, adenine-thymine and guanine-cytosine. We 
compare our results with available experimental data and 
previous calculations using methods from first principles. 

Subsection IIII Bl contains our calculations for the 
charge transfer hopping parameters between successive 
base-pairs (Sec. IHI B ip or between neighboring DNA 
bases fSec. IHI B 2\ . and provides estimates for the corre- 
sponding tight-binding parameters described previously 
in Section lll CI We present results for both electrons and 
holes and for all possible combinations of successive bases 
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TABLE I: For each DNA base (first row) are shown: the num- 
ber of atoms participating in tt bonds (second row) , the total 
number of contributed electrons (third row) , the ttHOMO 
energy ' (fourth row) and the ttLUMO energy ) 

(fifth row) of the isolated base, as well as the energy of the 



first TT-TT* transition E-j^- 



-,h (isol.) 



-E 



b {isol.) 
H 



(sixth row), 



G 



obtained in this work. We also show experimental values of 
the ionization energy IE (seventh row) and the first tt-tt* 
transition (eighth row), as well as existing theoretical predic- 
tions using methods from first principles for IE (ninth row) 
and -Ejr-ir* (tenth row). Finally, we list the ttHOMO and 
ttLUMO energies, E'^ (eleventh row) and E'l (twelfth row), 
respectively, of the bases when they are distorted within the 
base-pairs of B-DNA (see text), as obtained in our work. The 
quantities E^ in eleventh row and in twelfth row represent 
the parameters which appear in Eq. (|13p . All energies 



are given in eV. 



FIG. 1: (Color online) tt molecular structure of the 



bases: cytosine (C), thymine (T), guanine (G), and ^Mgteons 
(A). Full (empty) symbols correspond to occupied 
pied) molecular orbitals. 



or base-pairs. 



DNA base 



TT bonds 



isol.) 

ccu- 



E 
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TABLE IL For the two B-DNA base-pairs (first row) are 



shown: the ttHOMO energy E'^ (second row), the ttLUMO 



energy E'^ (third row) 
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HOMO and LUMO energies and wavetu nctioiis 
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Before starting our discussion about the HO.pgg^ipa 

LUMO of DNA bases, we show in Fig. [T] the if-ctcc- 
tronic structure of all bases, as obtained using thl§fc55[ 
method presented in Sec. Ill Al The geometric Hpjgti(jftJ3-DNA) 
ture of bases and the corresponding interatomic c iH^]^33-DNA) 
have been obtained from Ref. 44] . Occupied orbitals are 
shown by full symbols in Fig. [1] while unoccupied ones by 
empty symbols. The second and third row of Table U list 
-for each base of the corresponding column- the number 
of atoms participating in 7r-bonds and the total number 
of Pz electrons in the molecule, respectively. For exam- 
ple, adenine has twelve Pz electrons, which occupy, by 
pairs of opposite spin, the six lowest energy levels of A 
shown in Fig. [TJ 

Table U summarizes our results for the HOMO and 
LUMO energies of DNA bases. Relevant experimental 
values for the vertical ionization energy (which is equal to 
the absolute value of HOMO energy) and the first tt — tt* 
transition (i.e. the energy difference between LUMO and 
HOMO) are also shown, along with corresponding results 
of previous theoretical calculations using methods from 
first principles. When the bases form base-pairs within 
DNA, they are slightly deformed in comparison to their 
structure when isolated. The HOMO and LUMO ener- 
gies of the distorted bases in B-DNA (shown in the last 
two rows of Table |l| may differ from those of the iso- 
lated bases, and these values provide the relevant on-site 
energies E'^^j^ in the tight-binding Eq. (fT3l) . 
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and the corresponding first tt-tt* tran- 



sition energy Ett-tt* ~ E^ — Ejf (fourth row), as obtained 
in this work. We also list existing theoretical predictions us- 
ing methods from first principles, for the HOMO energy (fifth 
row) and the first tt-tt* transition (sixth row). The quantities 



E'^ in the second row and E"^ in the third row represent the 



parameters E'^^^^ 
given in eV. 
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FIG. 2: (Color online) Top: geometrical structure of ade- 
nine (right) and thymine (left), within an A-T base-pair. It is 
shown the standard numbering of atoms, as also used in this 
work. Cl' denote deoxyribose carbons and C5m thymine's 
methyl carbon, not participating in tt bonding. Middle: 
atomic occupation probabilities, |cip [cf. Eq. ([1])], for the 
HOMO and LUMO wavefunctions of isolated adenine (right) 
and thymine (left) bases. Bottom: atomic occupation proba- 
bilities, \Cif [cf. Eq. (O], for the HOMO and LUMO wave- 
functions of A-T base-pair. 



FIG. 3: (Color online) Top: geometrical structure of gua- 
nine (right) and cytosine (left), within a G-C base-pair. It is 
shown the standard numbering of atoms, as also used in this 
work. Cl' denote deoxyribose carbons, not participating in 
TT bonding. Middle: atomic occupation probabilities, |cip [cf. 
Eq. (m], for the HOMO and LUMO wavefunctions of isolated 
guanine (right) and cytosine (left) bases. Bottom: atomic oc- 
cupation probabilities, jCij^ [cf. Eq. dS])], for the HOMO and 
LUMO wavefunctions of G-C base-pair. 



In TableHIwe report results for the HOMO and LUMO 
energies (and the HOMO-LUMO energy gap) of the two 
B-DNA base-pairs, using the procedure desribed in Sec- 
tion lllB] Corresponding values obtained previously from 
first principles methods are also shown. 

Figure [2| (Figure [3]) presents for A, T, and A-T (for 



G, C, and G-C) the corresponding intramolecular occu- 
pation probabilities of holes or electrons, obtained from 
the calculated HOMO or LUMO wavefunctions, respec- 
tively. This is shown through the squared coefficients 
jcij^ (or jCij^ for base-pairs) in the expansion of Eq. ([T|) 
(or Eq. ([5]) for base-pairs) of the corresponding wavefunc- 
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tion. These squared coefScients represent the probabil- 
ities of finding an electron or hole (LUMO or HOMO, 
respectively) located at the corresponding atom of the 
base or the base-pair. The structures of the bases, pre- 
pared with the open-source software for chemical struc- 
tures Jmol [13], are shown on the top of each figure. The 
standard numbering of atoms has been used. 

A detailed discussion regarding the HOMO and LUMO 
wavefunctions is carried out below, separately for each 
isolated base and for the two B-DNA base-pairs. 



where the LUMO wavefunction is localized in atoms C8, 
C6, C2, N3, and N7. The only difference from our find- 
ings is the existence of a relatively higher probability in 
atom N7. In both results the lowest probabilities appear 
in the same atoms. Our findings are also in good ac- 
cordance with the results of Ref. 76] , where the LUMO 
wavefunction is localized in the same atoms (C8, 06, C2, 
N3). 

2. Thymine 



1. Adenine 

Our calculations show that the HOMO energy of ade- 
nine is —8.2 eV. This slightly underestimates the experi- 
mental value of the ionization energy (8.4 eV [S], 8.5 eV 
[isl . liol) by 0.2-0.3 eV. Our result is similar to that of 
Ref. [6J], where using ah initio calculations the HOMO 
energy was found —8.24 eV. Other ah initio calculations 
predicted the ionization energy between 8.3-8.6 eV [63l |. 
In good accordance with the measured values is also the 
ah initio prediction of 8.5 eV in Ref. |65|]. 

The calculated HOMO wavefunction (see Fig. [2l mid- 
dle right) is distributed all over adenine. In descending 
order, it has its main amplitude in atoms C5, C8, N3, 
N6, N7, C6, 02, 04, while the lowest probabilities exist 
in atoms Nl and N9. These results agree with those from 
Ref. [2l| , since in both cases the HOMO wavefunction is 
distributed all over adenine and the lowest probabilities 
are found at the same atoms. Our results are also in 
agreement with the outcome of density functional the- 
ory calculations [lH, where the tt HOMO wavefunction 
was found to be delocalized along adenine, with greater 
probabilities in atoms N7, N3, C6, 05, 08, and the lowest 
in atom N9, as well as with recent time-dependent den- 
sity functional theory predictions [t^ , where the HOMO 
wavefunction is distributed all over the molecule except 
from atom N9. 

We find the LUMO energy of adenine at —4.4 eV. Ac- 
cording to our calculations, the first tt-tt* transition en- 
ergy is E'tt-tt* = 3.8 eV. The deviation from the exper- 
imental value (4.5-4.8 eV [13, [H, H, [H, IM [sl) is 
0.7-1.0 eV. The i?,r-7r* value reported from several first 
principles methods is around 4.5-5.3 eV [sl, [tO, [ll, [tsI . 

i.e. the corresponding deviation from the exper- 
imental observations is 0-0.8 eV. It must be mentioned 
here that, although it is generally accepted for all DNA 
bases that the HOMO-LUMO energ y gap co rresponds to 
the first tt-tt* transition [H, HI, [ii, [sj, [S^l, in Ref. [77] 
the HOMO-LUMO gap of adenine is attributed to a n- 
TT* transition, while in Refs. [t^, [tI] the tt-tt* and tl-tt* 
transitions are predicted to be very close. 

Regarding the LUMO wavefunction of adenine, it can 
be seen from Fig. [2] (middle right) that the highest proba- 
bilities are found in atoms 06, 02, N3, and 08, while the 
remaining atoms are occupied by much smaller probabil- 
ities (< 5%). Our findings agree with those of Ref. [2l[, 



Thymine's HOMO energy is found —9.0 eV, which is 
in agreement with the experimental value of ionization 
energy, IEt =9.0-9.2 eV IH, gljp, [11] . The ah initio 
predictions of Refs. ^63] andTfUfall also within the 
range of experimental results, as they found the HOMO 
energy at —9.14 eV and —9.16 eV, respectively. However, 
other ah initio calculations [63| give the ionization energy 
between 9.4-9.7 eV, which is in worse agreement with the 
experimental data (a deviation of 0.2-0.7 eV). 

According to our calculations the HOMO wavefunction 
of thymine (cf. Fig. [2l middle left) is mainly localized at 
the atom 05 (with almost 50% probability) and at its 
nearby atoms 06 and Nl (with total probability around 
30%). The remaining atoms have probabilities less than 
10% each, with a negligible probability at the atom N3. 

We find that the LUMO energy of thymine is —4.8 eV. 
Looking at the first tt-tt* transition energy, our predic- 
tion of 4.2 eV underestimates the experimental results 
(4.6-4.7 eV [13, HH, [13, [H, [H, US Hi) by 0.4-0.5 eV. 
Calculated values from several ah initio methods 

are between 4.9-5.3 eV [H [H, [zl, [t^] , which overesti- 
mate the observed values by 0.2-0.7 eV. We mention that 
Refs. [13, [H III, [73 attribute the HOMO-LUMO gap of 
thymine to a ti-tt* transition. 

Fig. [2] (middle left) shows that the LUMO wavefunc- 
tion of thymine exhibits its main amplitudes in atoms 06 
and 04, with a total probability around 75%. All other 
atoms have less than 10% probability each, with vanish- 
ing amplitudes at the atoms 02 and 02. Our predictions 
agree with those from Ref. [tI], where the LUMO wave- 
function is mainly localized in atoms 04, 05, and 06. 



3. Guanine 

The HOMO energy of guanine is —8.2 eV according to 
our results. Our prediction for its ionization energy al- 
most coincides with the experimental observations (8.24 
eV [43] and 8.28 eV [48]). Previous ah initio calculations 
have found the ionization energy of guanine at 8-8.3 eV 
[H, 7.75 eV [li], and 8.21 eV [6|. These values deviate 
from the experimental ones by 0-0.5 eV. 

We observe in Fig. [3] (middle right) that the HOMO 
wavefunction of guanine is extended over a large part of 
the molecule, containing the atoms N2-C2-N3-C4-C5-N7- 
08 (with total probability around 90%). The atoms N3 
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and C5 present the highest probabihties (around 20% 
each), while about 15% probability exists in atom C8. 
A vanishing probability appears at the atom Nl and a 
very small one at N9. Our results agree with those of 
Ref. [81], where the main amplitude appears in the part 
N2-C2-N3-C4-C5 of the molecule and a lower probability 
also exists in atom C8. The main difference between the 
two calculations is the existence of a probability 18% at 
the atom 06 in Ref. [8l|, compared to 5% in our case. 
Similar to our findings are t he p redictions obtained from 
the ab initio method of Ref. [84! , where the HOMO wave- 
function is mainly localized at the same part of guanine, 
while there is also agreement with the results of Ref. [zl] , 
where the HOMO is predicted to be distributed all over 
the molecule apart from the atoms Nl, C6, and N9. 

We find that Guanine's LUMO energy is —4.4 eV. Re- 
garding the first TT-TT* transition, our prediction is 3.8 
eV, which underestimates the observed values (between 
4.3-4.5 eV |53,M,M,i9]) by 0.5-0.7 eV. The first tt-tt* 
excitation energy obtained from several first principles 
methods is between 4.4-5.3 eV (sl, [zl, [tI, [z3 > present- 
ing deviations of 0-0.8 eV from the experimental obser- 
vations. 

Looking at the LUMO wavefunction of guanine (Fig.[3l 
middle right), we see that the main amplitude appears in 
atom C6 with « 45% occupation probability. A signifi- 
cant probability also exists in atom C4 (more than 15%). 
All other atoms have probabilities less than 10%, with 
negligible amplitudes at C2, N2, C5, and N9. These re- 
sults differ from those of Ref. [84I , since there the LUMO 
wavefunction exhibits higher probabilities at the atoms 
C2, C4, Nl, C5, and N9, in descending order. On the 
contrary, the present findings are in reasonable agree- 
ment with those of Ref. UM ■ 



4- Cytosine 

The HOMO energy of cytosine is found —8.9 eV. Our 
prediction for the ionization energy coincides with the 
measured value of 8.94 eV [1^. Corresponding values 
calculated from ah initio methods are 8.87 eV [6J], 8.88 
eV [slj, and 9.0-9.4 eV The former predictions es- 
timate the experimental ionization energy equally well 
with our simple LCAO calculations, while the latter ones 
worse. 

It can be seen from Fig. [3] (middle left) that the HOMO 
wavefunction of cytosine is mainly localized at C5 (occu- 
pation probability « 45%), presenting also a significant 
amplitude at the neighboring atom C6 (probability larger 
than 20%). Atoms Nl, C2, 02, and N3 exhibit proba- 
bilities between 5% - 10%, while C4 and N4 have the 
smaller probabilities (below 2% each). In Ref. [8l| it was 
found that the atoms C5 and C6 accumulate only 25% 
probability, while the higher amplitudes appear in atoms 
02 and N3. 

We find the LUMO energy of cytosine equal to —4.4 
eV. The first tt-tt* transition energy according to our re- 



sults (4.5 eV) is compared well with the exp erimentall y 
observed values (4.5-4.7 eV [13, [H, [H, [EE [13, [H, [El, 
[60l . Frzj). Results from first principles methods predict 
values between 4.2-4.7 eV, i.e. some of them are within 
the range of experimental values [zl, [74 .FtsI. w hile others 
are slightly below the observed values 66l. l76l|. 

Figure [3] (middle left) shows that the LUMO wavefunc- 
tion of cytosine is mainly localized in atoms C6 and C4 
with occupation probabilities around 30% each. A prob- 
ability more than 10% appears in Nl, while C2, N3, and 
N4 have probabilities between 5% — 10%. These findings 
are in accordance with the results of Ref. [zl] , where the 
higher probabilities are found in atoms C4, C5, and C6. 



5. A-T base-pair 

As it has been already mentioned, the bases are slightly 
deformed within the base-pairs of B-DNA. This distorted 
form of the bases is appropriate when discussing base- 
pair properties (HOMO and LUMO energies or wave- 
functions), as opposed to the isolated form of the bases 
considered in the previous subsections. Therefore, be- 
fore going to the HOMO and LUMO discussion of base- 
pairs, we bricfiy mention the changes in bases' HOMO 
and LUMO due to the structural distortion within base- 
pairs. 

For adenine in the B-DNA conformation the HOMO 
energy is Ef:^ — —8.3 eV (differing by 0.1 eV in respect 
with the value of the isolated base, —8.2 eV, see Table|T|, 
while the LUMO energy remains unchanged {Ef^ — —4.4 
eV). Regarding the HOMO and LUMO wavefunctions 
of adenine the results are almost similar in both confor- 
mations. For example in the case of the HOMO wave- 
function, the highest difference of the occupation prob- 
abilities in respect with the results of the isolated base 
(Fig. [21 middle) is —1.2% in atom C2, while in all other 
atoms the differences are below 1%. In the case of LUMO 
wavefunction the differences are slightly higher, reaching 
in atom C8, -3.3% in atom C2, -2.4% in atom 
N3, and +2.1% in atom N7, while in the other atoms the 
changes are below 1.2%. 

For thymine the HOMO energy using the B-DNA 
structure is invariable in respect with that of the isolated 
base {EJj = —9.0 eV), while the LUMO energy changes 
slightly from —4.8 eV in the isolated form to — 4.9 
eV in the B-DNA conformation. The occupation proba- 
bilities for the HOMO (LUMO) wavefunction of thymine 
are almost similar in both cases and the differences are 
not larger than 0.5% (1.5%). 

According to our method, the HOMO energy of A-T 
base-pair is E^"^ = —8.3 eV. Comparing with ab initio 
results, this value differs only by 0.1 eV with the pre- 
diction of Ref. 64] (—8.19 eV), it is very close to the 
ones calculated in Refs. |zi] (-8.12 eV) and ^ (-8.06 
eV) , while there is a reasonable agreement with the value 
given in Ref. [zll (-7.8 eV). 

The LUMO energy of A-T base-pair is Ef^ = -4.9 
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eV. Therefore, the first tt-tt* excitation energy in our cal- 
culations is 3.4 eV. The prediction of i?,r-7r* obtained 
from first principles methods in Ref. [t^ is 6.39 eV. It 
must be mentioned that according to the latter study the 
HOMO-LUMO gap is not assigned to a tt-tt* transition, 
but to a TT-cr* transition [t^ . 

Figure[2] (bottom) shows that for the A-T base-pair the 
HOMO wavefunction is completely localized in adenine, 
while, on the contrary, the LUMO is completely localized 
in thymine. This is justified because of the significantly 
higher HOMO energy of A with respect to that of T 
(S^ - E]^ = 0.7 eV) and the significantly lower LUMO 
energy of T with respect to that of A — = —0.5 
eV), respectively, as compared with the corresponding 
overlap integrals of the bases' wavefunctions within the 
base-pair = 12 meV and |ti| = 9 meV, see Ta- 

ble |W1 below) . This also explains that E^'^ = E^ and 
E^^ = EJ^, as well as that the electron (hole) distribu- 
tion within the base-pair is identical with the correspond- 
ing distribution in thymine (adenine). The latter can be 
seen by a comparison of the corresponding wavefunctions 
in the middle and bottom panels of Fig. ^ taking also 
into account that adenine's HOMO and thymine's LUMO 
wavefunctions remain practically the same in the B-DNA 
and isolated conformations, as it is mentioned above. 



G-C base-pair 



The HOMO energy of guanine base in the base-pair 
conformation is E^ = —8.0 eV. There is a difference of 
0.2 eV compared to the corresponding result of the iso- 
lated base (-8.2 eV). Its LUMO energy in the B-DNA 
conformation {E^ = —4.5 eV) changes by 0.1 eV in com- 
parison to the value of —4.4 eV obtained for the isolated 
base. Concerning the HOMO wavefunction of guanine 
the highest differences of the occupation probabilities in 
respect with the results for the isolated base (shown in 
Fig. H middle) are -3.5% in atom N3 and +2A% in 
atom C8. In all other atoms the differences are below 
2%. For the LUMO state the biggest difference appears 
in atom C6, which is —4.2%, while in the other atoms 
the differences do not exceed 1.6%. 

Regarding cytosine, the HOMO and LUMO energies 
change slightly (by -1-0.1 eV), as compared to those of 
the isolated base. In particular, the predictions of cyto- 
sine's HOMO and LUMO energies within the base-pair 
conformation are E'^ = —8.8 eV and E^ = —4.3 eV, re- 
spectively. For the HOMO wavefunction, the differences 
of the occupation probabilities in the two base confor- 
mations are below 2.5% for all atoms. On the contrary, 
significant changes appear in the LUMO wavefunction of 
cytosine within the B-DNA conformation. The occupa- 
tion probabilities differ by -1-14.6% in atom C2, —6.4% 
in atoms C4 and C6, —4.1% in atom N3, and -1-3.2% in 
atom Nl, as compared to the isolated base LUMO, while 
in the other atoms the differences are below 2.5%. 

For the G-C base-pair, the obtained HOMO energy is 



E'^'^ = —8.0 eV. Several values calculated with methods 
from first principles have been reported in the literature: 
-7.68 eV [il, -7.51 eV [sil, -7.35 eV [zi], -7.2 eV 
[nl, and -6.27 eV 

The calculated LUMO energy here is = -4.5 eV. 
Concerning the first tf-tt* transition energy, our predic- 
tion is 3.5 eV. This result is in a reasonable accordance 
with the value obtained in Ref. [zl] (4.29 eV). 

However, these values differ from the one reported in 
Ref. [7§| (6.27 eV), where, similarly to the case of the 
A-T base-pair, the HOMO-LUMO gap is attributed to a 
TT-cr* transition. 

From Fig. [3] (bottom) we see that both HOMO and 
LUMO wavefunctions of a G-C base-pair are completely 
localized in guanine. The complete HOMO localization 
in guanine results again from the significantly higher 
HOMO energy of G in respect to that of C {E'^ - = 
0.8 eV), as compared with the hopping integral, |t/f | = 12 
meV in this case (see Table IIVI below) . Similarly, for 
the LUMO wavefunction of G-C, its complete localiza- 
tion in guanine is due to the lower LUMO energy of G 
with respect to that of C {Ef - e£ = -0.2 eV), as 
compared with the corresponding overlap integral of the 
bases' wavefunctions within the base-pair (t^ = 16 meV, 
see Table ITV)) . These results are in accordance with our 
findings that E^^ = E^ and E^^ = E^ . They also 
explain that the electron or hole distribution within the 
base-pair is the same with the corresponding distribu- 
tion in guanine, as it can be seen by a comparison of 
the corresponding wavefunctions in the middle and bot- 
tom panels of Fig. [3l Any small differences in the occu- 
pation probabilities of these plots reflect the changes in 
guanine's HOMO or LUMO wavefunctions between the 
B-DNA and isolated conformations, as it has been dis- 
cussed above. 

It is worth mentioning here that G-C's LUMO wave- 
function is the only HOMO or LUMO of both base-pairs 
which shows a qualitatively different behavior when the 
base distortions within B-DNA are taken into account, 
as opposed with the case of the isolated bases. In partic- 
ular, considering the structure of the isolated bases and 
calculating A-T's HOMO and LUMO or G-C's HOMO, 
one finds similar results with the ones presented in this 
and the previous subsection (for example, A-T's HOMO 
is locahzed in A and LUMO in T, while G-C's HOMO 
is localized in G). On the contrary, application of our 
method for the G-C LUMO, using the structure of the 
isolated bases, shows localization mainly in C (w 80%) 
but with a non- vanishing amplitude in G (« 20%). In 
this context note that the isolated G and C bases exhibit 
the same LUMO energy, -4.4 eV (see Table |l|. These 
results seem to be in accordance with the calculations of 
Ref. [3, although the nature of the HOMO and LUMO 
base-pair wavefunctions (whether it is of tt character or 
not) was not mentioned there. 
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B. Charge transfer hopping parameters 

The charge transfer hopping parameters tn and are 
obtained below, as described in Sec. |TT1 viz. through 
the HOMO and LUMO wavefunctions calculated in the 
previous subsection and Eqs. ^ or (fTTjl . for bases or 
base-pairs, respectively, with Vij given by Eq. 

1. Description at the base-pair level 

At first we present the hopping parameters t^^^ be- 
tween successive base-pairs, for the description at the 
base-pair level discussed in Sec. Ill C II (cf. parameters 

^ff/L^"^^^ in Eq. The notation YX is used here to 

denote two successive base-pairs, according to the follow- 
ing convention for the DNA strands orientation 

5' 3' 

^ ^compl 

3' 5' . (14) 

X, Xcompi, Y, Ycompi dcnotc DNA bases, where Xcompi 
and Ycompi are the complementary bases of X and Y, re- 
spectively. Therefore, the notation YX means that the 
bases Y and X of two successive base-pairs [Y-Ycompi 
and X- Xcompi) are located at the same strand in the di- 
rection 5' — 3'. X-Xcompi is one base-pair of the base-pair 
dimer and Y-Ycornpi is the other base-pair, separated and 
twisted by 3.14 A and 36°, respectively, relatively to the 
first base-pair. For example the notation AC denotes that 
the base-pair dimer consists of an adenine-thymine and 
a cytosine-guanine base-pair, where one strand contains 
A and C in the direction 5' — 3' and the complementary 
strand contains T and G in the direction 3' — 5'. 

Table IIIII summarizes the calculated hopping param- 
eters, according to our method, for all possible 
combinations of successive base-pairs, as well as corre- 
sponding values obtained from other theoretical studies. 
Due to the symmetry between base-pair dimers YX and 
Xcom,piYcompi, the number of different hopping parame- 
ters is reduced from sixteen to ten. In Table IIIII base- 
pair dimers exhibiting the same transfer parameters are 
listed together in the first column. Our predictions for 
the HOMO, ijf, and LUMO, t^l , hopping parameters 
are reported in the second and sixth column, respec- 
tively. The third column contains the ah initio results of 
Voityuk et al. |83(] for the magnitudes of HOMO trans- 
fer parameters. Results from Endres et al. [sj] for a few 
combinations of successive base-pairs, using the density- 
functional-based electronic structure program SIESTA, 
are presented for and t^l in the fourth and seventh 
columns, respectively. The same authors show similar 
results in Fig. 4 of Ref. The corresponding values 
-approximately extracted from that figure for zero twist 
angle- are shown in the fifth and eighth columns. As it 



has been already mentioned, in Ref. the transfer pa- 
rameters have been obtained using the same theoretical 
bedrock adopted also in our work -i.e. the Slater-Koster 
semi-empirical method and Eqs. ^ and (fTT|) - but with 
slightly different expressions for V^p^ and Vppa- 

TABLE III: Base-pair transfer integrals for all possible combi- 
nations of successive base-pairs, as shown in the first column 
according to the notation described in the text, in the direc- 
tion 5' — 3'. Hole hopping parameters t^^ obtained from our 
calculations (through HOMO wavefunctions), from Ref. [83| , 
from Ref. [sj, and from Fig. 4 of Ref. Q, are presented 
in the second, third, fourth, and fifth column, respectively. 
Electron hopping parameters tjf obtained from our calcula- 
tions (through LUMO wavefunctions), from Ref. [sj, and 
from Fig. 4 of Ref. 0), are shown in the sixth, seventh, and 
eighth column, respectively. These quantities represent the 
parameters t^y^'^*^' which appear in Eq. (|10|) . AH hopping 
integrals i^/^ are given in meV. 



Base-pair 
sequence 


.bp 
''H 


I4l 
[83] 


.bp 
''H 

[84] 


.bp 
''H 

f4] 


.bp 
''L 


.bp 
''L 

184] 


.bp 

f4] 


A A, TT 


-8 


26 


-70 


^ -25 


-29 


105 


^ 35 


AT 


20 


50 






0.5 






AG, CT 


-5 


122 


-71 


« -50 


3 


112 


« 35 


AC, GT 


2 


27 






32 






TA 


47 


55 






2 






TG, CA 


-4 


26 






17 






TC, GA 


-79 


25 


-187 


f« -160 


-1 


47 


^ 35 


GG, CC 


-62 


93 


-141 


« -140 


20 


53 


^ 35 


GC 


1 


78 






-10 






CG 


-44 


22 






-8 







We see from Table IIIII that according to our calcula- 
tions (second column), the higher in magnitude hopping 
parameters for holes appear in TC, GA, GG, CC, TA, 
and CG dimers, where the corresponding transfer matrix 
elements are around 50 meV or larger. The lower magni- 
tudes of hole hopping parameters (1-2 meV) occur in GC, 
AC, and GT dimers. According to Ref. [8^, the higher 
values are found in AG and CT dimers. In four out of 
sixteen dimers (GG, CC, CG and TA) the absolute val- 
ues of our predictions for and those of Ref. Q differ 
by no more than a factor of two, while in other five cases 
(GC, AC, GT, AG, and CT) our results are significantly 
smaller by more than one order of magnitude. However, 
it should be noted that more sophisticated theoretical de- 
scriptions [s^ have shown that the approximation used 
in Ref. 83] in general overestimates the transfer integrals. 
Endres et al. in Refs. [s^l and Q provide values for 
a few cases of successive base-pairs, which are larger in 
magnitude than both our results and those of Ref. [s^ , 
with the exception of AG and CT where their values are 
in between of our results and those of Ref. . 

Regarding electron hopping parameters (sixth col- 
umn), we find the highest magnitudes (around 30 meV) 
in AC, GT, A A, and TT dimers, while the lowest one is 
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in AT (less than 1 meV). The most recent t^^ values of 
Endres et al. in Ref. [1] are closer in magnitude to our 
predictions, as compared with their former estimates in 
Ref. ^8^. For the majority of base-pair dimers, the hole 
transfer parameters, as calculated in our work, arc larger 
in magnitude than the corresponding ones for electrons, 
indicating probably that hole transport is more favorable 
than electron transport along DNA. This is certainly true 
for the higher amplitude hopping parameters (larger than 
40 meV). 



2. Description at the single-base level 

At this level of description a charge carrier is allowed to 
hop not only between bases belonging to adjacent base- 
pairs, but also from one DNA base to its complementary 
base within the same base-pair [see the second term in 
the right-hand-side of Eqs. p3)) ]. In Table |TV] we display 
these interbase intra-base-pair hopping parameters. Our 
calculations show that the magnitudes of hole and elec- 
tron transfer integrals within G-C and A-T base-pairs are 
around 10 meV. For holes we find equal hopping param- 
eters within G-C or A-T base-pairs, while for electrons 
intra-base-pair hopping is more favorable in a G-C base- 
pair. First principles results from Ref. [s^ show larger 
magnitudes of electron hoppings than our predictions (al- 
most four times larger), while for holes the magnitudes 
are larger than ours in A-T and smaller in G-C base-pairs 
(in both cases the corresponding results differ by a little 
more than a factor of two). 



TABLE IV: Intra-base-pair hopping parameters of the two B- 
DNA base-pairs (first column), for holes {t\j, second and third 
column) and electrons (tj,, fourth and fifth column). Our 
predictions for holes (electrons) are shown in second (fourth) 
column, while results from Ref. [s^ are shown in third (fifth) 
column, respectively. These quantities represent the param- 



eters t 



i)(A,l;A,2) 
H/L 



and t 



6(A,2;A,1) 
H/L 



which appear in Eq. (|13p . All 



hopping parameters i^/^ are given in meV 



Base-pair 


T'H 












[86] 




[86] 


A-T 


-12 


26 


-9 


34 


G-C 


-12 


5 


16 


63 



Now we turn to interbase transfer integrals between 
bases of adjacent base-pairs. First we examine in- 
trastrand parameters, i.e. hoppings between successive 
bases of the same strand [see the third and fourth terms 
in the right-hand-side of Eqs. ((13])] . A similar convention 
in the notation is used here as previously; Y X implies 
that the DNA bases Y and X are given in the direc- 
tion 5' — 3' of the strand. For example, the notation 
TG denotes the strand orientation 5'— T— G— 3'. The in- 
trastrand interbase hopping parameters as obtained from 
our calculations are presented in Table |V] (in second 
column for holes and eighth for electrons). Regarding 



hole transfer, the higher in magnitude hopping matrix 
elements (larger than 100 meV) are found in TT and 
CT. Other relatively large magnitudes appear in TC, 
GC, GA, GT, AT, AC, CC, and GG (|4| > 60 meV). 
The lower absolute values of hole hopping parameters 
correspond to successive CG, AG, CA, and AA bases 
{\t\j\ < 10 meV). Looking at the electron hopping pa- 
rameters, the highest magnitudes appear in CT (larger 
than 60 meV) and CC, GC (|<^| > 40 meV), while the 
lower ones in AG, AC, AT, and TA {\t\\ < 10 meV). In 
most of the cases electron hopping parameters are smaller 
in magnitude than hole transfer parameters between suc- 
cessive bases, similarly with the case of hopping between 
successive base-pairs. 

Other theoretical predictions for intrastrand transfer 
integrals have been presented in Refs . [H, [H, [13, [H, Hlj , 
using methods from first principles. The corresponding 
values are shown in Table |Vl Because these methods are 
usually more reliable for occupied orbitals, there exist 
many works computing hole hopping parameters In 
general, our results are in the range of values obtained 
by these works. The ab initio method used in Ref. ^ 
(similar to that of Ref. [l^), has been shown to overesti- 
mate the transfer parameters [ssj. Indeed, in most of the 
cases the values provided in Ref. [sj are larger than the 
other values of Table [V] Regarding electron hopping 
parameters, first principles results of for a few cases of 



successive bases have been presented in Ref. 
are in a good agreement with our calculations 



which 



Finally we present interstrand transfer parameters de- 
scribing "diagonal" interbase hoppings between diago- 
nally located bases of adjacent base-pairs which belong 
to opposite strands [see the fifth and sixth terms in 
the right-hand-side of Eqs. (US])]- Here it is neces- 
sary to distinguish two sets of diagonal interstrand hop- 
pings, depending on the strand orientation. For the con- 
formation of two successive base-pairs shown in (jl4p . 
the one set of diagonal interstrand transfer parameters 
refers to hoppings in the direction 3' — 3' (i.e. the path 
3' — Ycompi ■ ■ ■ X — 3') and the other set refers to hoppings 
in the direction 5' — 5' (i.e. the path 5' — F • • • Xcompi — ^')- 



Table IVII shows results for interstrand transfer inte- 
grals in the direction 3' — 3' for all possible combinations 
of bases, while Tabic IVlTl presents corresponding results 
in the direction 5' — 5'. Our calculations yield relatively 
small magnitudes for these matrix elements. In general 
the obtained values are of the order of meV or tenths of 
meV, for both electrons and holes, apart from the cases 
of 3'-^ • • ■ A-3' and 3'-G- ■ • G-3' for holes (exhibiting 
between 40-50 meV) and 3' - A ■•• A - 3' for elec- 
trons (with around 30 mcV). Theoretical predictions 
for hole interstrand hopping parameters t''jj have been 
also presented for all combinations of bases in Ref. . 
Other results from fir st p rinciples methods obtained for 
a few cases in Refs. [83| (for holes) and [H^ (both for 
electrons and holes), show usually larger magnitudes in 
comparison to our calculations. 
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TABLE V: Intrastrand transfer parameters between two suc- 
cessive DNA bases (first column) in the direction 5' — 3' of 
the strand. Hole hopping parameters tj/ obtained from our 
calculations (through HOMO wavefunctions) , and from Refs. 
[stI ]. [ssI ]. [sa ]. [sgl ]. and liB], are presented in the second, 
third, fourth, fifth, sixth, and seventh column, respectively. 
Electron hopping parameters obtained from our calcula- 
tions (through LUMO wavefunctions) and from Ref. [80], are 
shown in the eighth and ninth column, respectively. These 
quantities represent the parameters ^^^j^'>-^='=i,») -^j^j^jj appgaj. 

in Eq. p3p . All hopping parameters tjj^/j^ are given in meV. 



Base 




141 


141 












sequence 




[87] 


[85] 


[88] 


[89] 


[86] 




[86] 


AA 


-8 


30 


4 


-4 


8 


21 


16 


25 


AT 


68 


105 




-63 


28 




7 




AG 


-5 


49 


44 


-10 


37 




1 




AC 


68 


61 




42 


30 




-3 




TA 


26 


86 




-31 


64 




-7 




TT 


-117 


158 




72 


93 


-98 


-30 


-23 


TG 


28 


85 


61 


18 


70 




-17 




TC 


-86 


76 




-28 


52 




22 




GA 


-79 


89 


36 


-77 


52 




30 




GT 


73 


137 


81 


141 


49 




-32 




GG 


-62 


84 


51 


53 


61 


-114 


20 


20 


GC 


80 


110 




-114 


57 




43 




CA 


5 


29 




-2 


5 




-12 




CT 


-107 


100 




-55 


33 




63 




CG 


-1 


42 




9 


31 




15 




CC 


-66 


41 




22 


26 


-21 


-47 


-60 



IV. CONCLUSIONS 

We have systematically studied all the tight-binding 
parameters which are necessary for the description of 
charge transfer afong DNA. The tt electronic structure 
of the four DNA bases (adenine, thymine, cytosine, and 
guanine) has been calculated by using the LCAO method, 
employed with a novel parametrization. In addition, we 
have presented the HOMO and LUMO of the two B-DNA 
base-pairs (adenine-thymine and guanine-cytosine) , us- 
ing a similar process of linear combination of molecular 
orbitals. Taking into account the slight deformation of 
bases within the base-pairs of B-DNA (compared to the 
isolated bases) we find that (a) for the A-T base-pair the 
HOMO resides in A and the LUMO in T, while (b) for 
the G-C base-pair both HOMO and LUMO reside in G. 

Our theoretical approach predicts the tt HOMO energy 
of DNA bases with a deviation smaller than 0.3 eV in 
comparison with the experimental values. In particular, 
for the bases of guanine and cytosine our results coincide 
with the experimental ones, while for thymine and ade- 
nine the difference ranges from 0-0.2 eV and 0.2-0.3 eV, 
respectively. Regarding the first tt-tt* transition energy 
the deviations are larger: 0-0.2 eV for cytosine, 0.4-0.5 



TABLE VI: 3' — 3' interstrand hopping parameters between 
two diagonally located DNA bases of adjacent base-pairs (first 
column). Hole hopping parameters t|f obtained from our 
calculations (through HOMO wavefunctions) and from Refs. 
[13 1 [13 1 presented in the second, third, and fourth col- 
umn, respectively. Electron hopping parameters t'j^ obtained 
from our calculations (through LUMO wavefunctions) and 
from Ref. [13 , are shown in the fifth and sixth column, respec- 
tively. These quantities represent the parameters f^^j^''^*^'-'' 
(i 7^ j) which appear in Eq. (|13|l . All hopping parameters 
are given in meV. 



Base 




4 


141 






sequence 




[88] 


[86] 




[86] 


AA 


48 


49 




29 




AT, TA 


-3 


-7 


9 


3 


-13 


AG, GA 


-3 


-11 




-6 




AC, CA 


-5 


17 




-3 




TT 


0.5 


6 




0.2 




TG, GT 


5 


-14 




2 




TC, CT 


0.5 


4 




-0.2 




GG 


-44 


-32 




-5 




GC , CG 


4 


22 


48 


-4 


-15 


cc 


1 


10 




0.3 





TABLE VII: 5' — 5' interstrand hopping parameters between 
two diagonally located DNA bases of adjacent base-pairs (first 
column). Hole hopping parameters t\j obtained from our 
calculations (through HOMO wavefunctions) and from Refs. 
[13 1 WA - [13 1 presented in the second, third, fourth, and 
fifth column, respectively. Electron hopping parameters 
obtained from our calculations (through LUMO wavefunc- 
tions) and from Ref. [s^ are shown in the sixth and seventh 
column, respectively. These quantities represent the parame- 
ters t^^y^'^^^'^'' {i 7^ j) which appear in Eq. (|13|) . All hopping 
parameters t\[^j^ axe given in meV. 



Base 


Iff 




141 








sequence 




[88] 


[87] 


[86] 




[86] 


AA 


2 


31 


31-35 




6 




AT, TA 


9 


7 


16-20 


-11 


2 


-10 


AG, GA 


4 


-13 


19-24 




3 




AC, CA 


5 


-1 






-2 




TT 


4 


1 


3-4 




2 




TG, GT 


5 


-9 






3 




TC, CT 


2 


0.3 






-2 




GG 


3 


12 


19 




-2 




GC , CG 


4 


2 




-8 


-3 


-12 


cc 


1 


1 


1 




2 





eV for thymine, 0.5-0.7 eV for guanine, and 0.7-1.0 eV 
for adenine. Our results for the HOMO and LUMO en- 
ergies and wavefunctions of bases and also of base-pairs 
have been compared with other theoretical calculations, 
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using methods from first principles. 

We provide estimates for the complete set of hopping 
parameters (both for electrons and holes) between succes- 
sive base-pairs and between neighboring bases (i.e., suc- 
cessive bases in the same strand, complementary bases 
within a base-pair, and adjacent base-pairs' diagonally 
located bases in opposite strands) in B-DNA, including 
all possible combinations of them. Our predictions for 
the transfer parameters arc compared with other theo- 
retical estimates, when available, and in most cases there 
is an agreement in the order of magnitude of the results. 
In general, the hopping parameters for hole transfer ob- 
tained through the HOMO wavefunctions, are higher in 
magnitude compared to the ones for electron transfer ob- 
tained through the LUMO wavefunctions. These theoret- 
ical calculations show that probably hole transport along 
DNA is more favorable than electron transport. 

The microscopic quantities calculated in this work, i.e. 



the HOMO and LUMO energies as well as the hopping 
matrix elements, provide all necessary parameters for 
a tight-binding phenomenological description of charge 
transfer along the DNA double helix, based either on 
bases' or base-pairs' tt molecular overlap. Taking advan- 
tage of such a description, at a mesoscopic level, the tem- 
poral and spatial evolution of electron or hole transport 
along DNA can be examined in experimentally relevant 
time and length scales. Furthermore, important measur- 
able quantities like hole/electron transmission coefficients 
and conductivities can be calculated for any DNA seg- 
ment (of whatever base sequence) under consideration, 
since there are well-established techniques to compute 
such a quantities from tight-binding models. 
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